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ABSTRACT 



We calculate radiative energy loss of heavy and light quarks in anisotropic medium (static) 
in first order opacity expansion. Such an anisotropy can result from the initial rapid lon- 
gitudinal expansion of the matter, created in relativistic heavy ion collisions. Significant 
dependency of the energy loss on the anisotropy parameter (£) and the direction of prop- 
agation of the partons with respect to the anisotropy axis is found. It is shown that the 
introduction of early time momentum-space anisotropy can enhance the fractional energy 
loss in the direction of the anisotropy, whereas it decreases when the parton propagates 
perpendicular to the direction of the anisotropy. 



1 Introduction 

One of the goals for the ongoing relativistic heavy ion collision experiments at the Relativistic 
Heavy Ion Collider (RHIC) and the upcoming experiments at CERN Large Hadron Collider 
(LHC) is to produce quark gluon plasma (QGP) and study its properties. According to 
the prediction of lattice quantum chromodynamics, QGP is expected to be formed when 
the temperature of nuclear matter is raised above its critical value, T c ~ 170 MeV, or 
equivalently the energy density of nuclear matter is raised above 1 GeV/fm 3 [1]. The 
possibility of QGP formation at RHIC experiment, with initial density of 5 GeV/ fm 3 initial 
density is supported by the observation of high px hadron suppression (jet-quenching) in the 
central Au-Au collisions compared to the binary-scaled hadron-hadron collisions [2]. Apart 
form jet- quenching, several possible probes have been studied in order to characterize the 
properties of QGP. 

However, many properties of QGP are still poorly understood. The most debated ques- 
tion is whether the matter formed in the relativistic heavy ion collisions is in thermal equi- 
librium or not. The measurement of elliptic flow parameter and its theoretical explanation 
suggest that the matter quickly comes into thermal equilibrium (with r thorm < 1 fm/c, where 
T therm is the time of thermalization) [3]. On the contrary, perturbative estimation suggests 
relatively slower thermalization of QGP [4]. However, recent hydrodynamical studies [5] have 
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shown that due to the poor knowledge of the initial conditions there is a sizable amount of 
uncertainty in the estimate of thermalization or isotropization time. It is suggested that 
(momentum) anisotropy driven plasma instabilities may speed up the process of isotropiza- 
tion [6], in that case one is allowed to use hydrodynamics for the evolution of the matter. 
However, instability-driven isotropization is not yet proved at RHIC and LHC energies. 

In absence of a theoretical proof favoring the rapid thermalization and the uncertainties 
in the hydro dynamical fits of experimental data, it is very hard to assume hydrodynamical 
behavior of the system from the very beginning. Therefore, it has been suggested to look 
for some observables which are sensitive to the early time after the collision. For example, 
jet-quenching vis-a-vis energy loss of partons could be an observable where the initial state 
momentum anisotropy can play important role. This is the issue that we address here. 

It is known that the energy loss of partons (also dubbed as 'jet-quenching') in QCD 
plasma can proceed in two ways: by two-body scattering and also via gluon radiation. 
These are known as collisional and radiative energy loss respectively. The phenomena of 
jet-quenching has been investigated by various authors [2] More recently, the non-photonic 
single electron data shows more suppression than expected which cannot be explained by 
radiative loss alone. A substantial amount of work has been done to look into this issue in 
recent times [2]. 

It is to be noted that the existing calculations on energy loss have been performed in 
isotropic QGP which is true immediately after its formation [7]. However, subsequent rapid 
expansion of the matter along the beam direction causes faster cooling in the longitudinal 
direction than in the transverse direction [4]. As a result, the system becomes anisotropic 
with (pl 2 ) << (pt 2 ) m the local rest frame. At some later time when the effect of parton 
interaction rate overcomes the plasma expansion rate, the system returns to the isotropic 
state again and remains isotropic for the rest of the period. Thus, during the early stage the 
plasma remains anisotropic and any calculation of energy loss should, in principle, include 
this aspect. The collisional energy loss in anisotropic media for heavy fermions has been 
calculated in Refs. [8, 9]. In these calculations it is found that the deviations from the 
isotropic results are of the order of 10% for £ = 1 and of the order of 20% for £ = 10. It is 
observed that the collisional energy loss varies with the angle of propagation by upto 50%. 

Recently in [10], the transport coefficient q has been calculated in anisotropic media, 
which in turn, affects the radiative energy loss. We here attempt to provide quantitative 
estimate of radiative energy loss by modifying the static scatterer model [11] appropriate for 
anistropic media. 

The other interesting aspect which, in recent years, has attracted considerable attention 
is the possibility of the growth of unstable modes in anisotropic plasma [12]. For example, in 
[13] the authors calculate q for a two-stream plasma and show that the momentum broaden- 
ing grows exponentially in time as the spontaneously growing fields exert an exponentially 
growing influence on the propagating parton. This momentum broadening of a fast parton 
which radiates gluons due to the scattering off the plasma constituents therefore controls 
the radiative energy loss [14]. In an evolving plasma this is an important component, which, 
however, is not included in the present manuscript. Therefore the results we report here can 
be considered to be something like a zeroth order approximation. 

The plan of the paper is the following. In section 2 we briefly mention how to calculate 
the 2-body potential in an anisotropic media along with the modified expression for the 
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fractional energy loss. Section 3 will be devoted to discuss the results. Finally, we conclude 
in section 4. 



2 Formalism 

In this section, we recapitulate the basic formalism of the radiative energy loss of a fast 
moving parton in an infinitely extended static isotropic plasma [15, 16, 11]. As in Ref. [11] 
we restrict ourselves to the radiative energy loss quarks at first order in opacity involving 
three diagrams as shown in Fig.(l) where we assume that an on-shell heavy quark produced in 
the remote past is propagating through an infinite QCD medium that consists of randomly 
distributed static scattering centers. In the original Gyulassy-Wang formalism [17] static 
interactions are modeled here as color-screened Yukawa potential originally developed for 
the isotropic QCD medium given by 

V n = V(q n )e l ^ x " 

= 27r5(q°)v(q n )e- i ^T an (R)®T an . (1) 

with v(q n ) = 4ira s /(qn + u 2 ), where \x is the Debye mass. x n is the location of the nth 
scattering centre, T (summed over a n ) denotes the colour matrices of the parton and the 
scattering centre. It is to be noted that the potential has been derived by using hard thermal 
loop (HTL) propagator in QGP medium. In a plasma with momentum anisotropy the two 
body interaction, as expected, becomes direction dependent. It has been observed that, on 
distance scale of the order of the inverse Debye mass, the attraction for the quarks aligned 
along the direction of anisotropy is stronger than for transverse alignment [18]. Therefore, 
the radiative energy loss, will also depend on the direction of momentum of the quarks 
emitting Bremstrahlung gluons. This necessitates the introduction of anisotropy dependent 
potential to estimate the radiative energy loss in a plasma having anisotropic momentum 
distribution. 
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Figure 1: Feynman diagrams contributing to the soft gluon radiation in a static medium to 
first order in opacity 



The heavy-quark potential in an anisotropic plasma has recently been calculated in [18] 
for which one starts with retarded gluon self-energy expressed as [20] 



UT{P) = g 2 J 



d 3 k u dfjk) 



v »pP * 

P ■ v + ie 



(2) 



We have adopted the following notation for four vectors: P M = (po,p) = (po,P,Pz), i- e. p 
(with an explicit vector superscript) describes a three-vector while p denotes the two- vector 
transverse to the z-direction. 

To include the local anisotropy in the plasma, one has to calculate the gluon polarization 
tensor incorporating anisotropic distribution functions of the medium. This subsequently 
can be used to construct HTL corrected gluon propagator which, in general, assumes very 
complicated from. Such an HTL propagator was first derived in [19] in time-axial gauge. 
Similar propagator has also been constructed in [18] to derive the heavy-quark potential in an 
anisotropic plasma, which, as we know, is given by the Fourier transform of the propagator 
in the static limit. 

The self-energy, apart from momentum P M , also depends on a fixed anisotropy vector 
n^{= (I, ft)) and Yi^ u can be cast in a suitable tensorial basis appropriate for anisotropic 
plasma in a co- variant gauge in the following way [18]: 

= aA^ + f3B^ + -fC^ u + 5D^ (3) 

where the basis tensors are constructed out of p M , n M and the 4-velocity of the heat bath w M . 
The detailed expressions for the quantities those appear in Eq.(3) can be found in Ref. [18]. 
The anisotropicity enters through the distribution function, 

f(p) = fU^ + ^V-nf) (4) 

where, the parameter £ is the degree of anisotropy parameter (—1 < £ < oo) and is given by 
£ = (p 2 (/(2)p 2 z )) — 1. It is to be noted that £ can also be related to the shear viscosity [21]. 

Since the self-energy is symmetric and transverse, all the components are not independent. 
After change of variables (p' — p*[l + £(p • n) 2 ] the spatial components can be written as 

11 _/i J 4tt V 1 + { + P-v + ie) (5) 

Now a, (3, 7 and 5 are determined by the following contractions: 

p*rrV = 

A il n l U ij p> = (f-(n-P) 2 )6 
A u n l U^A^ k n k = t^lll a + 1 ) 

TiU ij = 2a + /3 + 7 (6) 



where the expressions for a,/3,7 and 5 are given in Ref. [19]. 

After knowing the gluon HTL self-energy in anisotropic media the propagator can be 
calculated after some cumbersome algebra [10, 18]: 

4 2 \ 

+A G [(P 2 -a- i)^B^ + (co 2 - P)C^ + S^D^} - —P»P» (7) 
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where 



Aq 1 = (P 2 -a- 7 )(w 2 - f3) - 5 2 [P 2 - (n ■ Pf 



(8) 



Now the momentum space potential can be obtained from the static gluon propagator in the 
following way, 



where, 



m„ = 



m? = 



v(q,?*,0 = ^A uu ( W = 0,q,g,,O 

2 <f + m 2 t + m 2 
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g^tan- 1 ^)- 



(q 2 + m 2 + m 2 )((p + mj|) - m 2 



(9) 



= tan 1 



a/? 2 + £q 2 / 



/r 



(Vf + (1 + tan-i (v /£))(<f + Cq 2 ) + ^ tan- (-^L. 



2 



i + 2 4 

q 2 



£q 2 + 



. -i, fa ■ ^"W + 3£q 2 ) . -i/' 



*7J 2 £g*q|g1 
'4(£<f + <f) 3 / 2 



(10) 



with q = (q, g z ). For general anisotropy vector n we have, q = q — (g • n)n and q z = q- n. 
For, q z = 0, the potential in anisotropic media simplifies to 



q 2 + #(0/^ 2 



(11) 



+ 



tan" 1 



For small anisotropicity and q z = 
u(q,£ « 1) 



0, the two body interaction can be written as 
1 2 ^ 2 £ 



47ra s 
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q 2 + /i 2 3(q 2 + /i 2 ) 2 



(12) 



(13) 



Now in Fig. (1) the parton scatters with one of the colour centre with the momentum Q = 
(0, q, q z ) and radiates a gluon with momentum K = (uu, k, k z ). The method for calculating 
the amplitudes of the diagrams depicted in Fig. (1) is discussed in Refs. [15, 16] and we shall 
quote the main results only. The quark energy loss is calculated by folding the rate of gluon 
radiation (T(E)) with the gluon energy by assuming u> + go ~ In this approximation one 
finds, 



dE 
dL 



E r , dF 
— — / xdx— 
D R J dx 



(14) 
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Here Dr is defined as [t a , tj^cj^a] = Ci{G)CrDr where 02(G) = 3, Dr = 3 and [t a ,t c ] is 
a color commutator (see [11] for details), x is the longitudinal momentum fraction of the 
quark carried away by the emitted gluon. 
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Figure 2: Fractional energy loss for light quark when mean free is given is given by Eq.(19). 
where in anisotropic media we have, 



x 



dT 

dx 



C R a s L r d 2 k d 2 q 



71 



A 



7T 7T 



lQTT 2 ai 



k+ q 



(k + q) 2 + x 2 k 2 + X 



(15) 



In last expression, v(q, £) is the two body quark quark potential given by Eq.(9) and 



m 



2 x 2 + m 2 where m 2 = /i 2 /2 and m 2 = fi 2 /G. 



X = 

In the present case, we assume that the parton is propagating along the z-direction 
and the anisotropy vector n makes an angle 8 n with the z-axis, i.e. n = (sm9 n ,0,cos9 n ). 
Thus n describes the direction of propagation of the parton with respect to the anisotropy 



axis. In such case we replace q and q z in Eq. (9) by q — > J q 2 — q 2 sin 2 6 n cos 2 (j) and q z — > 
|q| cos^sin^n, where q = (|q| coscp, |q| sin0). 

For arbitrary £ the radiative energy loss can be written as 



AE 



C R a s Ljji 2 



7T Z 



A 



dxd 2 c{ 



+ 



l^(q ; g Z) 0l 2 r 1 

167r 2 a 2 \- 2 
q 2 - k 2 m + x 



+ X 2 J^ + 2cf(x-k 2 J + (k 2 m + x) 2 



+ 



q 2 + 2 X 



q- 



1 + ^ 



ln (" 



x 



: q 2 + 3 X ) + V / l + f(q 2 + x) 



X 



U2 



3x) + Ji + S vq 4 + 2 q 2 (A - k 2 J + (k 2 m + X ) 



(16) 
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Figure 3: (Color online) Same as Fig. 2 for £ = 1 (left panel) and for £ = 5 (right panel). 



In the above expression, A denotes the average mean free path of the quark given by 

1 1 1 

A Xg X q 



(17) 



which in this case would be £ dependent. In the last expression X g and X q correspond to the 
contributions coming from q-g and q-q scatterings. 
Explicitly with Eq.(ll) we have, 



X7 1 



C R C 2 {i)pi 



q 2 + i?(0/i 



2\2 ' 



where Cr = 4/3, 6*2(2) is the cashimir for ^-dimensional representation and C^i) = (N^ — 
1)/(2N C ) for quark and 62(2) = iV c for gluon scatterers. dA = N% — 1 is the dimensionality 
of the adjoint representation and pi is the density of the scatterers. Using pi = pf°/ y'l + £ 
we obtain 

1 _ 18a s TC(3) 1 1 + N F /6 

X 71*^1 + 1 R{i)l + N F /A 1 } 

where Np is the number of flavors. For £ — > Eq.(19) reduces to well-known results [11] 

18a s TC(3) 1 + N F /6 



1 

A 



1 + N F /A 



(20) 



It is evident that the changes from the isotropic media appear here as the co-efficient 
(-R(O) °f the Debye mass and the coefficient + of the number density. In the limit 

£ — > we recover all the previously known results as may be checked from Ref. [11]. 
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Figure 4: (Color online) Same as Fig. 2 for charm (left panel) and bottom quarks (right 
panel) with £ = 1 and T = 500 MeV. 



3 Results 

For the quantitative estimates of the fractional energy loss in an anisotropic media, first we 
consider a plasma at a temperature T = 250 MeV with effective number of degrees of freedom 
Nf =< 2.5 with the strong coupling constant a s =0.3 and L=5 fm. We also note that the mean 
free path of the propagating parton depends on the anisotropy parameter £ (see Eq.(19)). 
The fractional energy loss for non-zero £ (£=0.5) for light flavour is shown in Fig. (2). As is 
evident from Eq.(16), the energy loss in anistropic media depends on the angle of propagation 
of the fast partons with respect to the anisotropy axis (n). This is also illustrated in Fig. (2). 
It is observed that for non-zero value of the anisotropy parameter (£), the fractional energy 
loss increases in the direction parallel to the anisotropy axis. However, away from the 
anisotropy axis, the fractional energy loss decreases as the quark-quark potential is stronger 
in the anisotropy direction. 

For higher value of the anisotropy parameter £ the results are shown in Fig. (3). It is seen 
that the fractional energy loss increases with £ in the anisotropy direction. For £ = 1 and 
£ = 0.5, the fractional energy loss increases marginally for 9 n = ir/6 and it becomes larger 
for £ = 5 for the same value of 9 n . However, in the perpendicular direction the fractional 
energy loss decreases substantially. It is to be noted that for small anisotropy the 
results are almost similar to the case when the mean free path is independent 
of the anisotropy parameter. However, for larger values of £ the result changes 
reasonably as can be verified by calculating A from Eq.(19) for larger anisotropy 
(see Fig. (2)). 

For the heavy quarks, i.e. for charm and bottom, the results are shown in Fig. (4) for 
£ = 1. Similar to light quarks, we find enhancement in the anisotropy direction as well as 
for 9 n = 7r/6. However, for 6 n = tt/2 the energy loss (fractional) decreases for the reasons 
mentioned earlier. 
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4 Summary 



In this work, we have calculated the fractional energy loss due to gluon radiation in an infinite 
size anisotropic media treating the scatterer as providing a screened coulomb-like potential. 
We have seen that the potential gets modified in anisotropic media. It is observed that the 
fractional energy loss depends on the direction of propagation of the fast partons with respect 
to the anisotropy axis as well as on the anisotropy parameter £. An enhancement is seen 
in the direction parallel to the anisotropy direction ft, where as in the transverse direction 
it reduces due to weaker quark-quark interaction. It is also observed that for higher values 
of £, the fractional energy loss increases for a given direction with respect to the anisotropy 
axis. We also note that due to the dependency of the mean free path on the anisotropy 
parameter, the energy loss increases as £ increases. 

We do not include the recoil of the scatterer in this work. However, this condition can 
be relaxed by incorporating the recoil corrections which plays an important role as shown in 
Ref. [11]. This will be included in future publication. Furthermore, the finite size effect to 
the radiative energy loss in anisotropic media would also be interesting to study. 

The present calculation can be extended to include the effect of the growth of unstable 
modes to obtain results valid in a more realistic scenario as mentioned in the introduction. 
Inclusion of such effects might modify the quantitative estimate of nuclear modification factor 
at RHIC and LHC energies. 
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